Full counting statistics of transport through two-channel Coulomb blockade systems 
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A mesoscopic Coulomb blockade system with two identical transport channels is studied in terms 
of full counting statistics. It is found that the average current cannot distinguish the quantum 
constructive interference from the classical non-interference, but the shot noise and skewness are 
more sensitive to the nature of quantum mechanical interference and can fulfill that task. The 
interesting super-Poisson shot noise is found and is demonstrated as a consequence of constructive 
interference, which induces an effective system with fast-and-slow transport channels. Dephasing 
effects on the counting statistics are carried out to display the continuous transition from quantum 
interfering to non- interfering transports. 
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Rather than average current, the current fluctuations 
in mesoscopic transport can sometimes provide deep in- 
sight into the nature of transport mechanisms A 
fascinating theoretical approach, known as full counting 
statistics (FCS) theory [1, H, 

can simultaneously yield 
all the statistical cumulants of the number of transferred 
charges (i.e., all zero-frequency current-correlation func- 
tions). Experimentally, the real-time counting statistics 
has been carried out in transport through quantum dots 
[H, representing a crucial achievement of being able to 
count individual electron tunnel events. 

For charge transport at very low transmission, the 
uncorrelated transmission events are Poisson processes. 
However, the Fermi-Dirac statistics together with the 
possible many-body Coulomb interaction usually en- 
hances correlation among the transport electrons, and 
thereby results in sub-Poisson noise ||. It is thus of 
interest to examine mechanisms that can lead to super- 
Poisson-noise behavior, since the current fluctuations can 
be used in reverse to gain insight into the nature of 
unusual transport mechanisms. The mechanisms pro- 
posed so far for the super-Poisson noise include such 
as double electron charge transfer by Andreev reflection 
in normal-superconductor (NS) junction @, 0], multiple 
electron charge transfer by multiple Andreev reflections 
in SNS system [1, [t| [l(| QjJ , dynamical channel blockade 
0) EE EH , dynamical spin blockade [l5[ , bistability [lj| , 
cotunnelin g [171 fl8j , electron-phonon interaction in shut- 
tle system [19j, and decoherence in mesoscopic coherent 
population trapping system [20(. 

In this work we consider a relatively simple system, say, 
electronic transport through a Coulomb blockade system 
with two identical transport channels, which can be re- 
alized experimentally by transport through two adjacent 
levels in a single quantum dot (QD) [2l|, or through two 
QDs in parallel [12, HH . This type of setup itself is of 



particular interest, since it is an analogue of the opti- 
cal double-slit interferometer. In this context the un- 
derlying quantum interference and phase accumulations 
through QDs have been the subjects of intensive stud- 
ies [24], [H, . Our present study will be placed at the 
level of FCS, from which a number of interesting effects 
of quantum interference on current fluctuations will be 
revealed. 



Model Description 

The transport through a two-channel Coulomb block- 
ade system is governed by the Hamiltonian 



H — Hd + H Leads + Ht, 

Hd = Eid\d\ + E 2 d\d 2 + Unin 2 , 



(la) 



Hheads = '^2(£LkC Lk C Lk + £RkC' Rk C R k), (lb) 
fe 

H T = 5^[njtd}c ifc + n jR d]c Rk + H.c.]. (lc) 

jk 

Here c Lk Rk (cLk,n.k) and dl(dj) are the electron creation 
(annihilation) operators, for the electrode reservoirs and 
central dot states, respectively. The two channels are 
characterized by states with energy levels E\ and £2- 
Couplings of these two dot states to the electrodes are de- 
scribed by CljL(R), or = 2Trg LjR \Q jL j R \ 2 , for latter 
use. Here Ql,r are the density of states (DOS) of the elec- 
tron reservoirs. To manifest maximally the quantum in- 
terference effect, we shall focus on two identical transmis- 
sion paths. This can be accomplished by assuming equal 
and energy independent coupling strengths of the two dot 
states with the left and right electrodes, i.e., = 



n 



2L(R)\ 



Q 



L{R): 



and T l(R) 



V 2 

L L(R) 



L(R)- 
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address the quantum interference between transmissions 
through the two channels, the relative phase difference 
is significant. Physically, the phase difference contains 
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the phase accumulation of spatial motion from the elec- 
trode to dot, particularly in the presence of magnetic 
vector potential (i.e., the Aharanov-Bohm effect), as 
well as the phase changes associated with transmission 
through quantum dots [25|, [2(1 . Nevertheless, in this 
work we would like to adopt a phenomenological way to 
account for all these phase accumulations, by choosing 
= ^2Li an d rj — Hir/^r. Here rj can be regarded 
as a relative phase parameter. Note that the alternative 
gauge, say, fliR = fl^R and r\ — ili^/r^L, does not af- 
fect the final results [21(. In this paper, we also assumed 
77 = ±1, i.e., only complete constructive and destructive 
interference are considered. 

In the above Hamiltonian we omitted the spin indices, 
thus did not explicitly write out the on-site Coulomb in- 
teraction terms, and only left U n\ni to describe Coulomb 
interaction between electrons in the different dot states. 
In this work, unless explicit specification, our study will 
be restricted to the strong Coulomb blockade regime, 
which only allows for three available occupation states, 
i.e., |0), |1), and |2), corresponding to, respectively, 
empty dot state, and states with one electron on either 
E\ or E 2 . 



FCS Formulation 

Before going to the specific study of the above de- 
scribed system, we would like first to reformulate the FCS 
formalism based on the particle-number-resolved master 
equation approach [27|, |28i |29, |30| ■ As pointed out in Ref. 
|3U the pioneerin g w ork [2,3 and a few other approaches 
developed latter [id l32l. |33|] are largely restricted to ad- 
dressing the FCS of noninteracting electrons. While Ref. 
l3ll developed an elegant theory of FCS for mesoscopic 
systems in strong Coulomb blockade limit, however the 
system's internal quantum coherence did not enter it 
since the theory was constructed on the basis of classical 
stochastic processes. It is thus advantageous to formulate 
an approach of being able to account for both the internal 
quantum coherence and the many-body Coulomb inter- 
action on equal footing. Although such type of approach 
has been described and applied to coherent and interact- 
ing systems [H, 35] , for completeness we would like here 
to reformulate it, hopefully in a more transparent and 
unified way. 

To relate with our earlier work [2^, [3(|, we rename 
H s = H D , H B = H Le ads, and reexpress H' = H T as 

2 

H' = H T = ^2[d]Fj +R.c], (2a) 

F j = ^2 M jL C Lk + ^ SljRC Rk = f L j + f Rj . (2b) 

fc fc 

Regarding H' as perturbation, the second-order cumu- 
lant expansion leads to a formal master equation for the 



system's reduced density matrix [2{J [3(| : 

p(t) = -i£ P (t)- [ dT(c'(t)g(t,T)c'(r)gHt,T))p(t). 

Jo 

(3) 

Here the Liouvillian superoperators are defined as 
Cp = [H s ,p], C'p = [H',p], and Q{t,r){...) = 
G(t,r)(...)G ji (t,r), with G(t 7 r) the usual propagator as- 
sociated with system Hamiltonian Hg. The reduced den- 
sity matrix p(t) = TrB[pT(t)], and (...) = Tr b[---Pb] with 
Pb the density matrix of the electron reservoirs. 

The trace in Eq.([3]) is over all the electrode degrees 
of freedom. To achieve the FCS of current fluctuations, 
we would like to keep track of the records of electron 
numbers emitted from the source lead (ni) and arrived 
at the drain lead (n^)- We therefore classify the Hilbert 
space of the reservoirs as follows: B^ ni ' n2 ^ — B L ni ' ® 
B^ 2 \ The entire Hilbert space of electron reservoirs is 
thus decomposed as B — @ ni Jl2 B^ ni ' n2 \ 

With this classification the average over states in the 
entire bath Hilbert space in Eq. © is replaced with the 
average over states in the subspace £?( ni >™ 2 ), leading to a 
conditional master equation 

P {n ^\t) = -iLp^ n2 \t) - [ drTr B(ni ,„ 2) 

Jo 

[C'(t)Q(t,T)C'{T)G\t,T)p T (t)}. (4) 

Here, p( ni > n2 \t) = Tr B („ 1 ,„ 2) [p T (i)], is the reduced den- 
sity matrix of the central system conditioned by the elec- 
tron numbers emitted from the source lead (ni) and ar- 
rived at the drain lead (/12) until time t. 

To proceed, following Ref. [13, two physical consider- 
ations are further implemented: (i) Instead of the con- 
ventional Born approximation for the entire density ma- 
trix pr(t) — p(t)(£)pB, we propose the ansatz pr(t) — 
E ni> „ 2 P (ni ' n2 \t)® p { ^- n2 \ where p { ^ n2 \t) is the den- 
sity operator of the electron reservoirs associated with 
ni-electrons emitted from the source and ^-electrons 
entered the drain. The orthogonality between reservoirs 
states in different subspaces leads to the term selection 
from the entire density operator pt- (ii) Due to the closed 
nature of the transport circuit, the extra electrons en- 
tered the drain will flow back into the source via the ex- 
ternal circuit. Also the rapid relaxation processes in the 
reservoirs will quickly bring the reservoirs to local ther- 
mal equilibrium characterized by the chemical potentials. 
As a consequence, after the state selection procedure, the 
electron reservoirs density matrices p( ni '" 2 ) should be re- 
placed by . 

Further use of the Markov-Redfield approximation 
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leads Eq.(|l]) to an explicit form: 

7j ip (ni.»2) = J^idj^ P (ni ' n2) + P (ni, " 2) 4 +) 4 

- A^V" 1 '™ 2-1 ^] - d}p^' n2+ ^A^)] + H.c. . 

(5b) 

Here A%> = E, A$ = 

and = J2 a =L r ^aj ■ The spectral functions 

are defined in terms of the Fourier transform of 
the reservoir correlation functions, i.e., cj^j(±£) = 

f°° dtC^f} (t)e ±lCt . The reservoir correlators read 

{fUWaiir)) = C%](t - r), and (/ m (t)it,(r)) = 

C^](t-r). Here (■ ■ •) stands for Tr B [- ■ • p { °\ with 
the usual meaning of thermal average. Obviously, 
(Fjmir)) = Cfit-r) = E^L, R C^(t-r), 

and (Fi(t)f}(r)) = C[j\t-r) = E a=L . R C$ (t - r). 
For the sake of brevity, the explicit expressions of the 
reservoir correlation functions, the corresponding spec- 
tral functions, and A^ are ignored here, and are pre- 
sented alternatively in Appendix A. 

At this stage, it is worth making a few remarks as 
follows: (i) The above particle- number-resolved master 
equation is applicable to finite temperatures, which is an 
extension of Gurvitz's approach 27]. (ii) The second- 
order cumulant expansion of H' restricts the applica- 
bility to the regime of sequential tunneling. However, 
generalization to higher order expansion of H' [H, [37j 
and self-consistent corrections [38j are possible. The cor- 
responding FCS version is an interesting subject for fu- 
ture work, (hi) The above (m, ri2)-resolved master equa- 
tion generalizes the result in Ref. 30, for electron count- 
ing from at one junction to at two junctions. Further 
generalization to multi-terminal setup is straightforward, 
following precisely the same treatment, (iv) The con- 
nection of the particle-number-resolved density matrix 
with the distribution function of FCS is obvious, i.e., 
P(m,n2,t) = Tr[p(t)( ni > n2 '], where the trace is over the 
central system states. From this distribution function, 
all orders of cumulants of transmission electrons can be 
calculated. 

In practice, instead of obtaining the distribution func- 
tion from the solution of the particle-number-resolved 
master equation, a more efficient method is the cumu- 
lant generating function (CGF) technique. In the fol- 
lowing study, we only consider single counting statistics. 
That is, we only keep n>2, after making summation over 
n\. Multiple counting statistics in multi-terminal setup 
follows the same technique. 

Mathematically, the CGF is introduced as 

e -F( x ) =Y J P{n,t)e inx . (6) 



Here \ corresponds to the so-called counting field. 
Based on the CGF, the fe t h cumulant reads Ck — 
— {—id x ) k F(x)\x=o- For instance, the first two cumu- 
lants give rise to the mean value of the transmitted elec- 
tron numbers C\ = n, and the variance C2 = n 2 — n 2 ; 
the third cumulant (skewness), C3 = (n — n) 3 , charac- 
terizes the asymmetry of the distribution. Here, (■■■) = 
Eti("')^(' 1 '0- Moreover, the cumulants are straight- 
forwardly related to the transport characteristics, e.g., 
the average current by / = eC\/t, and the zero- frequency 
shot noise by S — 2e 2 C2/t. The Fano factor is defined 
as F — C2/C1, which represents the amplitude of cur- 
rent fluctuations, with F > 1 indicating a super-Poisson 
fluctuation, and F < 1 a sub-Poisson process. 

The generating function can be calculated as follows. 
Define S(x,t) = E„P (n) 0)e mx - Obviously, e~ F W = 
Tr[S(x, *)] ■ Let us reexpress the particle-number-resolved 
master equation formally as 

= V l) + Cp {n+l) + Dp^-V, (7) 

then S(x, t) satisfies 

S = AS + e~ ix CS + e lx DS = C X S. (8) 

The formal solution reads S(x,t) — e Cxt S(x,Q)- In the 
low frequency limit, the counting time is much longer 
than the time of tunneling through the system. One 
can prove [H EH, HI , that F( X ) = -\i(x)t, where 
Ai(x) is the eigenvalue of C x , and satisfies the condition 
Ai(x) lx-o- 0. 

FCS Analysis 

We now turn to the specific system under study. In 
the strong Coulomb blockade regime and at zero temper- 
ature, the matrix element form of the particle-number- 
resolved master equation reads 

p^ = ^Tj^+Tnip^+P^) 

+ v T R (p[ n r 1) +P 2 n r 1) ), (9a) 

.(n 2 ) r< ("2) t-i ("2) *-R 1 (n 2 ) . (na)\/nu\ 

Pn = t lPoq - - 7-j—{p\ 2 + p 21 ; )(9b) 

•("2) t-i ("2) T-i ("2) R / ("2) 1 ("2) \ fn \ 

P22 = T LPoq - T R p 22 - V—(Pl2 + P21 )>( 9c ) 

p^ = iSep^ + T^-Tnp^ 

-V^(p { n 2) + P^ 2) )- (9d) 

Here we have summed n\ and remained only 712, indicat- 
ing the mere study of FCS of the electrons entered the 
drain reservoir. For clarity, we have denoted the level 
spacing by 6e = E2 — El, and choose the reference of 
zero energy such that E2 — Se/2, and El = —Se/2. In 
the derivation of Eqs. ©, we assumed that E\^ are in- 
side the window of bias voltage, i.e.,/iz, > £1,2 > PR, and 
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U is infinite. In Ref. l2l| the same master equation was 
derived by a wavefunction approach. 

Performing a discrete Fourier transformation 



to Eqs.©, we obtain 
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-77^ 


i<5e — Ti? 
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r L 





-iSe - T R J 



(10) 

According to the definition of the cumulants we can ex- 
press Ai(x) as 



fc=i 



k\ ' 



(11) 



Then insert the above expansion into |Ai(x)/ — C x \ = 0, 
and expand this determinant in series of [ix] ■ Since 
the value of i\ is arbitrary we can obtain Ck by setting 
the coefficients of (ix) h equal to zero and solving them 
sequentially [2(j. Analytic expressions of the first two 
cumulants are accordingly obtained as 



2T L T R 6e 2 



(12a) 



C 2 



2r L (Se2-(r 1 -l)r R )+r R 6e 2 ' 
5e 4 T 2 R + 4r|[fe 4 + 2jeV| + (77 - 1^T|] 

[2r i (5e 2 -(7 7 -l)r|)+r fl <5e 2 ] 3 
x2r L r fl< 5e 2 . (12b) 

while the higher order cumulants can be instead carried 
out numerically, to avoid their lengthy expressions. 

In Fig. 1 the first three cumulants of transport cur- 
rent are displayed. It is of interest to note that in the 
Coulomb-blockade regime a super-Poisson noise is devel- 
oped by the constructive interference between the two 
paths (i.e. 77 = +1). This is clearly shown by the solid 
curve in Fig. 1(a). With the increase of the coupling 
asymmetry (i.e. a = Tr/Fl), the super-Poisson feature 
will be more evident. In contrast, for destructive interfer- 
ence (77 = —1), the current fluctuation is sub-Poissonian, 
as plotted by the dashed curve in Fig. 1(a). 

Another intriguing finding is that the super- and sub- 
Poisson characteristics are associated with different be- 
haviors of the skewness C3/C1 , as shown in Fig. 1 (b) . For 
destructive interference, the skewness is approximately 
zero, meanwhile for constructive interference, transition 
of the skewness from (small) positive to (large) negative 
values takes place, by increasing the coupling asymmetry 
(a). As is well known for photon counting statistics in 
quantum optics, the skewness (both its magnitude and 
sign) provides further information for the counting statis- 
tics, beyond the second order cumulant. As a compari- 
son, the results of non-interacting system are plotted in 
Fig.l (c) and (d), where neither the super-Poisson noise 
nor the negative skewness is found. 

To understand better the above super-Poisson be- 
havior, below we present an analysis in terms of fast- 
and-slow transport channels. Let us introduce an 




FIG. 1: First three cumulants of the zero- frequency current 
fluctuations. The solid and dashed curves display the results 
for constructive (rj — 1) and destructive (77 = —1) interfer- 
ences, respectively. The variable a = Tr/Tl characterizes 
the asymmetry of dot-state couplings to the left and right 
electrodes, (a) and (b) show the Fano factor and skewness 
for strong Coulomb blockade system, while (c) and (d) for 
system allowing double occupancy, i.e., U = 0.0, as a com- 
parison. Inset of (b) plots the average current / = eCi/t, 
with the convention e = 1. In the calculation Fl = l.OSe was 
assumed. 



alternative representation for the double dot states, 
with the corresponding electronic operators f± — 
(flLdx + n R d 2 )/y/n 2 L + fi£, and h = {^Ld 2 - 



toRdi)/y/fi 

-S, 



fl R , as well as the state energies Ex/ 2 = 



^y. In such representation the entire Hamiltonian is 
reexpressed as 

H = E x flfx + E 2 flh + 7 (/t/ 2 + /t/j) 
+ y^[On / / 1 t ai,fe + Oih/Jcir/s + H.c] 



(13) 



+ yfozLfzOLk + & 2 nfia Rk + H.c] 



Note that the strong Coulomb blockade, rather than be- 
ing explicitly described in this Hamiltonian, is reflected 
alternatively by the single occupation of the two dot 
states. For the sake of brevity, explicit expressions for 
the coupling between the new dot states, and their cou- 
plings to the electrodes are presented in Appendix B. 

In the new states representation, the formation of the 
fast-and-slow transport channels is demonstrated in Fig. 
2: (i) for constructive interference (77 = +1), increasing 
a = T R /T L can lead to F iL > T 2L , and T 1R » T 2R ; (ii) 
for destructive interference (77 = —1), however, increasing 
a = T R /T L leads to T 1L > T 2L , but T 2R »T 1R . Obvi- 
ously, in the former case, an effective fast-and-slow trans- 
port channels are evaluated, but can not in the second 
case. Such fast-and-slow transport channels will lead to 
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FIG. 2: Effective couplings of the dot states to the electrodes 
in the transformed state representation. It is found that ef- 
fective fast-and-slow transport channels are evaluated for the 
constructive interference (rj = 1), which is contrasted remark- 
ably with the destructive interference (rj = —1). Inset: effec- 
tive coupling between the transformed dot states. Parameters 
a = T R /F L , and F L = 6e = 1.0. 



bunching behavior, and result in the super-Poisson noise 

0, E sum. 

Decoherence Effect 

To illustrate the decoherence effect, we consider de- 
phasing between the two interfering paths, which is de- 
scribed by including the matrix element form of the fol- 
lowing Lindblad-type terms 

into Eq.(9), where the jump operators L^i = \/T2\l)(l\, 
and 1/02 = \/Td|2)(2|. The effects of decoherence on the 
first three cumulants are shown in Fig. 3, where the solid 
and dashed curves correspond to the constructive (r) = 1) 
and destructive (77 = —1) interferences, respectively. In- 
terestingly, for the constructively interfering transport 
(j] = 1), dephasing does not influence the transport cur- 
rent, see the solid line in Fig. 3(a). This is in remark- 
able contrast with the result of double-slit optical inter- 
ference, where the constructively interfering intensity is 
four times of the intensity of the individual path (slit), 
while the non-interfering intensity is simply two times 
of the single-slit intensity. This essential difference is 
originated from the multiple forward-and-backward scat- 
tering between the dot states and the electrodes in the 
case of electron transport. However, also for 77 = 1, 
the second and third cumulants (i.e. C2 and C3) sen- 
sitively depend on the dephasing strength. In particular, 
dephasing would cause a transition from super-Poisson 
to sub-Poisson processes, meanwhile the skewness (C3) 
changes from negative value to zero. The present analy- 
sis clearly shows that the super-Poisson current fluctua- 
tion is a consequence of the constructive interference. For 
destructively interfering transport (77 = —1), the almost 
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FIG. 3: Dephasing effects on the first three cumulants, which 
show a continuous transition from quantum interfering to 
classical non-interfering transports. The solid and dashed 
curves display results for the constructive (7/ = 1) and de- 
structive (77 = —1) interferences, respectively. Parameters 
a = T R /Ti = 3.0, and F L = Se = T. 



vanished transport current will be restored by dephas- 
ing, whereas the shot noise and skewness approximately 
do not change with dephasing. 



Conclusion 

To summarize, we have presented a FCS study for 
transport through a mesoscopic Coulomb blockade sys- 
tem with two identical transport channels. The FCS 
analysis showed that the shot noise and skewness are 
more sensitive than the average current to the under- 
lying quantum interference. In particular, the average 
current cannot distinguish the quantum constructive in- 
terference from the classical non-interference, while the 
shot noise and skewness can fulfill that task. In the 
regime of quantum constructive interference, the inter- 
esting super-Poisson shot noise was found, and was un- 
derstood in terms of an effective fast-and-slow channel 
picture via state-representation transformation. Dephas- 
ing effect has also been carried out to display the continu- 
ous transition of the first three cumulants from quantum 
interfering to classical non-interfering transports. Exper- 
iments within current technology are capable of examin- 
ing the predictions of this work. 

Methodologically, the present work extended the or- 
thodox FCS approach in strong Coulomb-blockade sys- 
tems [3lj . The present approach can not only account 
for the many-body Coulomb interaction, but also handle 
the internal quantum coherence. This advantage allows 
for a wider range of applications. 
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APPENDIX A 

In this appendix we present the explicit expressions 
of the spectral functions in Eq. (5) in the main text. 
For noninteracting electrodes and under the wide-band 
approximation, the reservoir correlation functions simply 
read 



(Ala) 

C$(t-T) = 9R&R I dsRke^^n^isRk), 

(Alb) 

Cmi(t-T) = V9R^R J dsRke^^-^n^ieRk) 

= c£\(t-r) (Ale) 

Here, the indices i and j denote the dot-states |1) and |2). 
Since the two dot-states are almost degenerate in energy 
and are equally coupled to the electrodes, in the above 
we have assumed that VL ia Q* a = Q* a Qj a = £1^. Note 

that an exception is Q^f^ij = Q^R^-iR = *7^ij> where 
i] = ±1 is the relative phase factor. This is because 
we have attributed the phase difference to the coupling 
amplitudes with the right electrodes. The electron and 
hole occupation functions are introduced, respectively, as 
n { a\e ak ) = n a (e ak ), and n { a\e ak ) = 1 -n a (e ak ), with 
n a {e a k) the Fermi function. Fourier transformation of 
Eq. (Al) gives the spectral functions: 



(A2a) 
(A2b) 

(7^(w)=» 7 r fl n^ ) (Tw), (A2c) 



-(Ej + Urij)dj, the op 



where T a — 2irg a n^. 
Furthermore, using Cdj 

erators in Eq. (5), which is defined by A^ = 

Ei C%){+C)di, A£ = C^}{-C)di, are accordingly 
obtained as 



A 



(+) 

aj 



-(Ei + Un^di 



APPENDIX B 

Via the transformation of state-representation as de- 
scribed in the main text, the energy levels of the trans- 
formed dot-states and their effective coupling strength 
read 

^_6e Sl* R -ill 



2 n\ + n R ' 

- _ge Oj-Q| 
_ Ql^r 



(Bl) 



Simple algebra also gives rise to the effective coupling 
strengths of the dot-states with the electrodes: 

Q 1L = ,^ 1L =a (CIl + Qr), 



QlR 
^2R 



tMT&r 

&-2R _ 
&-2R 



(nn L + n R ), 
{n L - n R ), 

(fl L - T)il R ). 



(B2) 



More transparently, for T] = 1, the corresponding tunnel- 
ing rates read 



Til 

F 2 R 



r L (r L + r fl + 2VrZr^) 

r L + r R ' 

r R (r L + r R + 2yTTO 

Tl + Tr 



r L (r £ + r fl -2yrZr^) 

r L + r R ' 

rfl(r L + r^-2VrZr^) 



(B3) 



T L + r R 

and for rj = — 1, they are 

T l (T l + T r + 2VTlT^) 



IL 



Tir = 



(A3a) 



r fl (r L + r fl -2Vnr^) 

r L + r R : 

r L (r L + r fi -2VrZr^) 

T2L ^TT~r ' 

ri?(r L + r fi + 2vTOT) 



(B4) 



2R 



T L + r R 



.(-) 



Y J C { J[HE l + Un^]d l - (A3b) 



Here the index i simply means differing from i, i.e., i = 2 
if i = 1, and vice versa. 
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